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We discuss electric transport and orbital magnetism of multilayer graphenes in a weak-magnetic 
field using the matrix decomposition technique. At zero temperature, the minimum conductivity 
is given by that of the monolayer system multiplied by the layer number A*', independent of the 
interlayer hopping t. When the interlayer hopping satisfies the condition t ^ h/r with r being 
collision time of impurity scattering, [N/2] kinks and [N/2\ + 1 plateaux appear in the Fermi- 
energy (gate voltage) dependence of the conductivity and the Hall conductivity, respectively. These 
behaviors are interpreted as multiband effects. We also found that the Hall conductivity and the 
magnetic susceptibility take minimum value as a function of temperature, for certain value of the 
gate voltage. This behavior is explained by Fermi-energy dependence of these functions at zero 
temperature. 
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I. INTRODUCTION 



Experimental studies of graphene have revealed ex- 
otic transport properties such, as an anomalous quan- 
tum Hall effect and the finite universal conductiv- 
ity at zero energyiii^ These results are essentially ex- 
plained by the two-dimensional (2D) massless Dirac 
equation which describes the low-energy band struc- 
ture of graphene around the gapless point i^i^i^i^i^i^ 
Moreover, multilayer graphenes which consist of stacked 
few-layer systems also attract attention. For bilayer 
systems, quantum Hall effccti'^'iS and longitudinal 
conductivit y ^^'^^i^^i^^ have been studied. For systems 
with more than three layers, electronic structures are in- 
vestigated experimentallyi^'i^iii and theoreticallyj ^^'^^i^° 
One of the most interesting point of these multilayer 
systems would be variety of stacking structures. A 
graphene is usually produced by micromcchanical cleav- 
age of graphite, so that the stacking structure is consid- 
ered to be the Bernal type, since the natural graphite 
falls into this category. However, production of graphene 
with other stacking types may also be possible by re- 
cent epitaxial methods. The difference of band structure 
depending on the stacking types are discussed )^^i^° and 
stability of the stacking structures is also studied in terms 
of symmetry arguments 

On the other hand, the magnetic susceptibility of car- 
bon systems has been studied for long times, before 
the discovery of graphene. It is well known that three- 
dimensional graphite shows large diamagnetism, and this 
has been explained theoretically by McClure. He showed 
that the orbital diamagnetism appears with a delta func- 
tion peak at the zero energy point based on the 2D mass- 
less Dirac equationj ^^i^^ This argument can be essentially 
applied to monolayer graphcnesi^ and effects of impu- 



rity scattering2^i^ and of an energy gap21 are discussed. 
The orbital magnetism in multilayer systems was studied 
more than two decades ago, motivated by graphite inter- 
calation compounds. Especially, the magnetic suscep- 
tibility of bilayer and multilayer systems was discussed 
by Safrau)^ and by Saito and Kamimura^ respectively. 
Quite recently, Koshino and Ando calculated the sus- 
ceptibility using matrix decomposition techniquci^ They 
discussed that Hamiltonian of the Bernal stacking sys- 
tems can be block diagonalized into effective bilayer and 
monolayer Hamiltonians depending on parity of layer 
numbers. This is a powerful tool to investigate multi- 
layer systems. 



In this paper, we turn our attention to the electric con- 
ductivity and Hall conductivity of multilayer graphenes 
in a weak-magnetic field. We also consider how the differ- 
ences of stacking structures appear in the physical quan- 
tities. We use the matrix decomposition technique used 
by Koshino and Ando throughout this paper for the cal- 
culation of Bernal stacking systems. We also discuss the 
finite-temperature properties including magnetic suscep- 
tibility in these systems. 



The rest of paper is organized as follows. In Sec. [U 
we discuss the Hamiltonian of the multilayer systems. In 
Sees. |TTT] and |TVl the conductivity and the Hall conduc- 
tivity arc discussed, respectively. Kinks and plateaux ap- 
pearing in their Fermi-energy dependence are discussed. 
In Sec. |Vl we discuss the finite-temperature properties 
of the diamagnetic orbital susceptibility, and discuss a 
minimum value as a function of temperature. In Appen- 
dices, we discuss the decomposition of Hamiltonian, and 
present analytical forms of physical quantities. 
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(a) Bernal stacking 






(b) rhombohedral stacking 




FIG. 1: (Color online) Two stacking structures of multilayer graphene. 
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FIG. 2: (Color online) Band structure of multilayer systems (A^ = 1-5) with the Bernal stacking [Eq. ((T|, solid line] and with 
the rhombohedral stacking [Eq. ([2} , dashed line] . For the Bernal stacking, linear dispersion relations appear for odd layer cases. 
The dispersion relation of the effective Hamiltonian of the rhombohedral stacking [Eq. dotted line] is also shown. 



II. STACKING STRUCTURES 



given by 



As is well known, the band structure of graphene has 
two gapless points due to the hexagonal lattice, and the 
low-energy property of this system is described by the 2D 
massless Dirac equation. For the multilayer graphene, 
we consider network of the Dirac fermion systems con- 
nected by interlaycr hopping t. It is also known that 
there are mainly two stacking types for graphite: one is 
the Bernal (staggered) stacking where the layer sequence 
can be written 1212 ... , and the other is rhombohedral 
stacking 123123 .... (see Fig. [J) 



Hh = 



7r+ 
t 



7r+ t 

t 7r_ 



(1) 



The Hamiltonians of the Bernal stacking systems arc 



7r_ 
t 



where tt± = vtt± with v being the velocity of the Dirac 
equation for the monolayer system, and tt± = i itTj,. 
TT = p + eA /c is the momentum operator in a magnetic 
field W X A ~ {0,0, B). The band structures obtained 
as eigenvalues of this Hamiltonian are shown in Fig. [51 
Quite recently, Koshino and Ando proved that the effec- 
tive Hamiltonian of A^-layer Bernal stacking system be- 
comes, without loss of generality, isolated [A^/2] bilayer 
system {[x] is the maximum integer which docs not ex- 
ceed x) with effective interlaycr hoppings tl,t2, ■ ■ ■ , ^ [)v/2] ' 
and one monolayer system if N is odd.— This result 
can also be obtained by calculating determinant of the 
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Schrodingcr equation, (see Appendix [X| The values of 
effective hoppings are given in Table HI We will use this 
decomposed effective Hamiltonian throughout this paper. 
For the bilayer system, an energy gap A = t between 
two bands appears, since the energy spectra are given by 
£ = ±(v/i2 + {2vhkY ±t)/2. 

On the other hand, the rhombohedral stacking is de- 
scribed by 



t 



7r+ 

t 



7r+ 
t 



7r+ 



(2) 



The band structure of this system is complicated, as 
shown in Fig. [2l and decomposition such as the Bernal 
stacking is difficult. If we turn our attention to the two 
bands near the zero-energy point, the effective Hamil- 
tonian becomes the following 2x2 formiSiSl (see Ap- 
pendix E]) 



1 



N-l 



(w^+)" 
(w7r_)^ 



(3) 



This effective Hamiltonian is useful to discuss the 
quantum Hall effect^^ and the zero-energy longitudinal 
conductivity ; "'^^'^"^ but it is difficult to obtain physically 
relevant results for the Hall conductivity and the suscep- 
tibility in the weak-magnetic field treatment below (see 
Appendix [B|. Therefore, calculations in this paper is 
mainly devoted to Bernal stacking systems. 

The low-energy effective Hamiltonian ([3]), however, 
gives information of the reason why the Bernal stacking 
is more stable than the rhombohedral stacking: Since the 
dispersion relation is e = zt{vhk)^ /t^~^ , the density of 
states of the A^-layer rhombohedral system is given by 



DNie) = 



Vt 



2/N-l 



(4) 



where V is volume of the system. For > 3, Dm{s) di- 
verges at e = 0, where the density of states of the Bernal 
case is I?2(e) = const. Therefore, the rhombohedral sys- 
tems are considered to be unstable against external per- 
turbations. 



III. CONDUCTIVITY 

First, we consider the conductivity of the multilayer 
systems based on the linear response theory. The con- 
ductivity is given by the Kubo formula^ 



Re Uui, = lim 



Im n^^(0,cj + iO) 



(5) 
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TABLE I: Effective interlayer hopping integral t* (i = 
1, 2, . . . , [Af/2]) of the effective bilayer Hamiltonian of the 
Bernal stacking systems. 



where H^^(q',Lj) = n^^(q,w) - H^^(q,0). The polariza- 
tion function in the Matsubara form is given by 

^f^M.^^rrd = y j dTe"'"^(r,J,,(q,r)J,(O,0)), (6) 



n^i/(o,iz^,„) 



^^tr (^Tm^'+Ti^): 



(7) 



k n 



where /? = l/k^T and V are the inverse temperature 
and the volume of the system, respectively. The current 
operator is given by 

M-q, r)^-eJ2 + 9/2, r)7p*(fe - q/2, r), (8) 



where ^{k,T) is the Fourier component of the field opera- 
tor ^'t(r) = [■0]i^(r),-0]3^(r), V]i,(r),'0tj^(r), . . .], where 
Ai and Bi indicate two sublattices of the hexagonal lat- 
tice of zth layer. The matrix 7^ is defined by 



I on 



(9) 



In Eq. (I7|), the impurity-averaged temperature Green 
function is given by 

g = g{k,lUJn) = {lUJn + [fi + lSgn{uJn)r-nQ]/h)-\ (10) 

where w„ = (2?! + l)7r//3ft is the Matsubara frequency 
of fermions. For = Q[k^uJn + ^'m), i^m = 2-Kml(3h 
is the Matsubara frequency of bosons. Here, we have 
introduced the scattering rate phenomenologically as the 
quasiparticle self-energy F = — Im E^, neglecting the fre- 
quency and the momentum dependences. This parame- 
ter stems from the scattering of the impurity potential 
implicitly assumed in the present system, and is related 
to the mean free time of quasiparticles as F = h/2T^ 

The final analytic form of Uxx for a bilayer system is 
presented in Appendix IB] In Fig. [31 the zero-temperature 
conductivity (calculated per valley and per spin) a^x ver- 
sus the Fermi energy /i scaled by the scattering rate F is 
shown. Experimentally, is a tunable parameter by the 
gate voltage. At zero energy /i/F = 0, the conductivity 
takes the minimum value axx = Ne^/irh which is the 
N times of the minimum conductivity of the monolayer 
systemj^i^i^ii For this finite value cr,„in = /irh, there are 
controversial argumentsy^i but all the analytic results do 
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not coincide with the experimental value e except 
for the numerical analysis of Nomura and MacDonald.— 
As the Fermi energy increases as ^l/T —t oo with small 
interlayer hopping t/T 1, the conductivity approaches 
to N times of the Drude-Zener form, 



1 + {uJ.Tf 



(11) 



Compared with the result of the linearized Boltzmann 
equation for nonrelativistic electrons do = ne^r/m, 
where n is the electron density, the "electron mass" m is 
related as m = l/zj/w^ii This agreement with the Boltz- 
mann description is due to the constant F. A similar 
linear Fermi-encrgy dependence of the conductivity is re- 
produced by numerical calculation considering effects of 
screened Coulomb impurity scattering.— 

For large interlayer hopping i/F ^ 1 which means 
that the energy gap A ~ 0{t) satisfies the condition 
A ^ K/t, kinks appear in the Fermi-energy dependence. 
The number of kinks is increased as the layers are in- 
creased. These phenomena can be interpreted as a multi- 
band effect: When the interlayer hopping is large enough 
in the iV-layer system, N bands in the positive (nega- 
tive) energy region split into [{N + l)/2] gapless modes 
and other bands with different energy gaps. Therefore, 
the number of kinks in the positive or negative energy re- 
gion is N /2 for even layers, and [N — l)/2 for odd layers, 
reflecting the discontinuity of the density of states. 

In order to clarify this argument, we generalize Eq. (jlip 
to the multiband system for // > as 



o N 



(12) 



i=l 



where the velocity Vi and the Fermi wave number of the 
i-th band i are related as 



hkp^i, /-t — ei{kY,i). 



(13) 



In Fig. [31 the Fermi-energy dependence of the zero- 
temperature conductivity at ^/F = 6 for various layer 
numbers is shown. The results of scmiclassical anal- 
ysis show good agreement with those of the linear re- 
sponse theory, except for the small energy region where 
the quantum effect is essential. For more detailed semi- 
classical argument, we should consider the Boltzmann 
equation for Dirac type systemsi^i 



IV. HALL EFFECT 



Next, we consider the Hall conductivity Uxy in a weak- 
magnetic field. The expression of the Hall conductivity in 
terms of the Green function is obtained by the Luttinger- 
Kohn representation^ for the basic functions and the 
Fourier expansion of the vector potential A{r) = AqC^'^'^ . 
As the first-order perturbation of the current term of the 
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FIG. 3: (Color online) Fermi-energy (/^/F) dependence of 
the longitudinal conductivity of Bernal stacking systems with 
= 2, 3, 4, 5, scaled by the minimum conductivity of the 
monolayer system cTmin = /nh. In these results, [A'^/2] kinks 
appear for the positive (or negative) Fermi-energy regions. 



Hamiltonian — Aq ■ J{—q)/c, we have the following 



5 



40 
30 

-g^ 20 
10 





quantum / / 




semiclassical / 




^N=4 / / / 



















-15 -10 -5 5 10 15 

n/r 



of hall coefficient p^y = '~o-xy/[{o'xx)'^ + {o'xy)'^] at the 
surface of a graphite»2£ 

The most remarkable feature of the Hall conductivity 
is that plateaux appear when the interlayer hopping is 
large t/T ^ 1. The number of plateaux for the positive 
or negative energy regions is (TV + 2)/2 for even layers, 
and {N + l)/2 for odd layers, and axy is "quantized" by 
the unit of WcTCTq. This phenomenon is due to effect of 
the energy gap A ^ 0{t) satisfying the condition A 3> 
K/t. Actually, a plateau also appears in the result of 
the monolayer system with an energy gapi^ In order to 
clarify this behavior, we also generalize the Drude-Zener 
theory for axy in the multiband systems. Then the Hall 
conductivity for /i > is given by 



FIG. 4: (Color online) Comparison between the longitudi- 
nal conductivity obtained by the linear response theory (solid 
lines) and that by the semiclassical analysis (dashed lines) for 
one to five layer systems with t/T = 6. 
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X e"''"^(r, J^g, r) J„(-g, r') J.(0, 0)). 



(14) 



Then by q expansion of the temperature Green function 
with the relation d^^Q = Q"tfj,G, we obtain the polar- 
ization function in the linear order of the magnetic field 
B = i{q,Ay - qyAl) as 



H^t.(0,u/„) = 1 



e^B 1 
2ch V(3h 



E 

k.n 



X tr (7^(?+7^C;+7^C;7j^(7 - 7^g_^7y(7+7^g7^g 

+ lt,G+lxG+lyG+luG - lt,G+lyG+lxG+lvG 

+ itiG+iuGixGiyG - liiG+iuGiyGixG)- (15) 

In Fig. [5l we show the zero-temperature Hall conduc- 
tivity (calculated per valley and per spin) Uxy versus the 
Fermi energy (gate voltage) scaled by the scattering 
rate F. The Hall conductivity changes the sign depend- 
ing on the sign of p.. For small t/T, the Hall conductivity 
shows sharp change in the small energy region /i/F — > 0. 
Then its absolute value takes maximum, and approaches 
to the constant value (oc F~^), as the Fermi energy is 
increased as /i/F — > oo. This value is N times of the 
Drude-Zener like formula, 



N 



(17) 



where the velocity Vi and "electron mass" m are given by 
Eq. (fT3|) . In Fig. [6l the Fermi-energy dependence of the 
zero-temperature Hall conductivity axy at t/T = 6 for 
various layer numbers is shown. The result of the semi- 
classical analysis shows good agreement with those of the 
linear response theory in terms of the step structure of 
axy- For odd layers, the peak structure remains even 
when the plateaux appear, because of the contribution 
from a gapless monolayer mode. 

We also show in Fig.[7l the Hall conductivity for = 3 
rhombohedral stacking system calculated based on the 
Hamiltonian ([2]). Reflecting the difference of the band 
structure (see Fig. [5]) , the first plateau appears at axy — 
ztWcTCTo, while axy — -iiluJcTao for Bernal stacking. 



V. ORBITAL MAGNETISM AND 
FINITE-TEMPERATURE PROPERTIES 

The orbital magnetism of multilayer graphene with 
Bernal stacking at zero temperature has already been 
discussed by Koshino and Andoi^ Here, wc turn our at- 
tention to the finite-temperature properties of the mag- 
netic susceptibility. The general formula for the orbital 
magnetic susceptibility of Bloch electrons was derived by 
Fukuyama;^ A modified expression of this formula which 
is also applicable for the effective Hamiltonian of Eq. ([3]) 
(see Appendix [B]) is 



n k 



(18) 



1 + {uJctY 



(16) 



where lUc = \eB\/mc is the cyclotron frequency. The de- 
viation of the Hall conductivity from this classical value 
at the low-energy regions and change of the sign are often 
called "anomalous" Hall effect. Actually, this behavior is 
observed experimentally in the gate voltage dependence 



where 7± = (7^; =F Yyy)/2. This result is obtained by 
the Luttinger-Kohn base, the Fourier expansion A{r) = 
Aq(e"' '' — e~"^ '')/2i, the q expansion of the temperature 
Green function, and the second derivative of the thermo- 
dynamic potential x 

field is given as B 
with q 
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qxAi- 
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where A^ 
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qyA(=il/2){q+. 
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FIG. 5: (Color online) Fermi-energy {fJ./T) dependence of 
the Hall conductivity of Bernal stacking systems with A'^ = 
2, 3, 4, 5, scaled by the classical Hall conductivity of the mono- 
layer system WcTcto. In these results, [A'^/2]-|-l plateaux appear 
in large t/T cases. 



Figure [S] shows the susceptibility of the bilayer sys- 
tem at zero temperature, (see Appendix [Bl) As discussed 
by Safrau)^ the susceptibility diverges logarithmically at 
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FIG. 6; (Color online) Comparison between the Hall conduc- 
tivity obtained by the linear response theory (solid lines) and 
that by the semiclassical analysis (dashed lines) for one to five 
layer systems with t/T = 6. 
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FIG. 7: (Color online) Fermi-energy (/i/F) dependence of 
the Hall conductivity of a rhombohedral stacking system 
with TV = 3. In this case, the first plateau appears at 
a^y — izujcTao, refiecting the band structure. 



zero-energy point with F = 0, and paramagnetic regions 
appear. As t/T is increased, the Lorentzian-like curve 
becomes broad, and the diamagnetism becomes smaller 
(larger) at small (large) Fermi-energy region. This be- 
havior is similar to that of the monolayer system with 
an energy gap.— Therefore, as discussed by Koshino and 
Ando, odd layer systems with Bernal stacking tend to 
show large diamagnetism than the even layer systems 
due to the contribution form gapless mode of the effec- 
tive monolayer. 

In Fig. [9l we show the temperature dependence of the 
susceptibility of monolayer and bilayer systems. We find 
that there appears the minimum value at finite tempera- 
ture when the gate voltage satisfies the relation with some 
critical value \fi\ > ^c- This phenomenon is explained in 
the following way: For simplicity, let us consider the sus- 
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FIG. 8: (Color online) Magnetic susceptibility of the bilayer 
system at zero temperature as a function of the Fermi en- 
ergy jj, for (a) fixed interlayer hopping t = 1 and for (b) 
fixed scattering rate F = 1. These data are scaled by 
Xo = —Qe^v^ /tv^c^ . t dependence of the susceptibility is sim- 
ilar to that of gap dependence in the monolayer system [see 
Fig. 1(a) of Ref.[23. 



ceptibility in the monolayer system with finite F,— 

1 



X 



da;/(.T)Im 



{x + ir)2 ' 



(19) 



where f{x) = (e*^^^ + 1) ^ is the Fermi distribution 
function. At zero and high temperature limits, Eq. (jl9p 

2 2 p 

becomes x = --ij^-p-pr^ and x = 0, respectively. In 
the finite-temperature regions, Sommerfeld expansion of 
(Ha is 



where 



F 



+ — ikBTyF'{fi) + 0{kBTy, (20) 



F2 ' 6 



F{x) = Im 



{x + iF)2 ■ 



F'(x) = ^ 21 



Therefore, for F' {[i) < {\fi\/r > l/VS) a mini- 
mum value appears as a function of temperature. This 
means that the condition to appear a minimum value 



FIG. 9: (Color online) Temperature dependence of the mag- 
netic susceptibility x of (a) the monolayer system and (b) the 
bilayer system for several strength of the interlayer hopping 
at /i/r — 1.5. A minimum value appears as a function of 
temperature. 



is existence of an inflexion point in the susceptibility 
at zero temperature as a function of the Fermi-energy 
oc F{x)dx. The same argument can also be applied 
to multilayer systems and other physical quantities, and 
it turns out that the Hall conductivity may have mini- 
mum value for small interlayer hopping or odd layer cases 
(see Fig. [T0|) . where the peak structure of the Hall con- 
ductivity is clear. Similar argument for a minimum value 
of physical quantity at finite temperature is recently done 
for the magnetization of quantum spin chains in a mag- 
netic field4i For the longitudinal conductivity, there ap- 
pears no minimum (see Fig. Ilip . 



VI. SUMMARY AND DISCUSSION 

In summary, we have studied the electric transport and 
the orbital magnetism of multilayer graphene in weak- 
magnetic field. We have found that kinks appear in 
the gate voltage dependence of the conductivity, and 
plateaux in the Hall conductivity. These phenomena are 
explained as multiband effects which become clear when 
the energy gap and the collision time satisfy the condition 
A ^ H/t. We have also considered finite-temperature 
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FIG. 10: (Color online) Temperature dependence of the Hall 
conductivity a^y of (a) the monolayer system and (b) the bi- 
layer system for several strength of the interlayer hopping at 
fj,/T = 3. A minimum value appears as a function of temper- 
ature. 



FIG. 11: (Color online) Temperature dependence of the con- 
ductivity axx of (a) the monolayer system and (b) the bi- 
layer system for several strength of the interlayer hopping at 
fj,/r = 3. They behave monotonically as a function of tem- 
perature. 



properties of this system, and found that a minimum 
vahie appears in the magnetic susccptibihty and the Hall 
conductivity, as functions of temperature. These phe- 
nomena are explained by the existence of an inflexion 
point in the zero-temperature Fermi-energy dependence. 

In this paper, we have tuned our attention mainly 
to the Bernal stacking systems, except for the result 
of Fig. [71 For rhombohedral staking systems, the band 
structure is quite different from that of the Bernal stack- 
ing, especially the gapless linear dispersions for odd TV > 
3 do not appear. Therefore, it is expected that difference 
of the properties of the physical quantities, such as the 
magnetic susceptibility, is not so drastic depending on the 
parity of the number of layers. It is also reported that the 
trigonal wrapping effect causes essential difference in the 
longitudinal conductivity^^ This would be also an inter- 
esting future problem. We have also discussed the semi- 
classical analysis of transport properties using the results 
of Boltzmann equation for nonrelativistic electrons. It is 
expected, however, that Boltzmann equation for Dirac 
type systems is needed for more detailed description in 
regions close to the zero-energy point. Moreover, it is 
also desirable to extend the present analysis to finite- 
magnetic field regions where the Landau quantization is 



essential. 

Note added in proof: After the submission of this pa- 
per, we become aware of a preprint H. Min and A.H. 
MacDonald, arXiv:0711.4333 where the present matrix 
decomposition technique and the 2x2 effective Hamil- 
tonian are used to discuss the quantum Hall effect in 
multilayer graphenes with general stacking structures. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN 

In order to calculate physical quantities of multilayer 
systems, it is useful to reduce the original Hamiltonian 
with 2N X 2N matrix form into some effective Hamil- 
tonian which has fewer matrix elements. McCann and 
Fal'ko reduced the 4x4 matrix Hamiltonian into the 2x2 
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form to discuss the quantum Hall effect of the bilayer sys- 
tem, which describes the two bands near the zero-energy 
point.— In order to derive such effective Hamiltonian, we 
should find Ticff which have the same eigenvalues of the 
original Hamiltonian 7i: 



dct(e-7i) =dct(e-Hc 



0. 



(Al) 



For the rhombohcdral iV-laycr system the 2x2 
effective Hamiltonian is obtained in the following way: 
First, we change the order of the matrix elements 
fmmAi,Bi,A2,B2,... to . . . S3, ^2, ^i, ^i, ^2, A3, . . . , 
where Ai and indicate two sublatticcs of the hexago- 
nal lattice of ith layer. Then, we have 



t 



t 



fc_ t 
k+ t 
k_ t 
k+ t 



(A2) 

where k± = i iky is eigenvalue of the momentum 
operator 7r_|-. We have set v = h = 1 for simplicity. Next, 
we calculate det(e — H) approximately for e t, then 
we obtain 



det{e-n) 



k 



2N 



0. 



(A3) 



One of effective Hamiltonians which satisfies Eqs. (|Aip 
and l|A3[) can be chosen as Eq. ([3]). 

On the other hand, the effective Hamiltonian for the 
Bcrnal stacking graphenes is discussed by Koshino and 
Ando, quite rcccntlyi^ According to their result, A^-layer 
Bernal stacking system can be described by isolated [A^/2] 
bilayer systems with some effective interlayer hopping, 
and one monolayer system if N is odd. This is exact 
mapping of Eq. ([T]) without using any approximation. 
Here, we derive the same result by the argument based 
on the determinants. After the same reordering of the 
original Hamiltonian ([1]) as the rhombohcdral stacking, 
we have 



k. 



k+ 

k^ t 
k+ t t 

k_ t t 

t '■■ 



(A4) 



Then we obtain the following recursion relation for the 
determinant of the A^-layer system A]y = det(e — H) as 



= _ 1^2^ 



A 



N 



(A5) 
(A6) 

(A7) 



Note that recursion relation such a closed form cannot be 
obtained for the rhombohcdral stacking. Equation (|A7p 
can be solved exactly, and factorized in the following way: 

Al =X, 
A2=X^ - t^e^, 
A3 ==X{X^ - 2fe'^), 



13 = 

A4 = 



V5-l \ ^2^2 



A5 = X{X^ - t^e^){X^ - 3t^e^), 



Thus, the determinant for the iV-layer system can be de- 
composed into those of the bilayers [X^ — (t*)^e^] with 
the effective hopping t* (Table |T| and that of one mono- 
layer X for odd A^. This means that the original Hamil- 
tonian ([1]) can be block diagonalizcd into subsystems. 
Koshino and Ando have obtained the effective hopping 
as t* = 2tsin[rmT/2{N + 1)] with m being an appropri- 
ate integer, by considering diagonalization of the matrix 
() A4p .— This factorization of the Hamiltonian is analo- 
gous to the A'^-leg ladder in quantum spin systems where 
a flat dispersion appears for odd N,^ 



APPENDIX B: ANALYTIC RESULTS 

Here, we present anaclitic forms of physical quantities 
discussed in this paper. Wc also briefly discuss the re- 
sults based on the low-energy effective Hamiltonian of 
the rhombohcdral stacking ^ given by 2 x 2 matrix. 

The longitudinal conductivity is given by the following 
form, 



dxYf'{x)]AN{x), 



(Bl) 



where /(x) — (e'^'^ + 1) ^ . For the monolayer case 
— 1, we have 



Ai(a;) 



(B2) 



For the bilayer system N = 2 with 4x4 Hamiltonian, 
the result is 
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-^2(2;) = -—-5 <^ 77 - Taii"^ — } + — — J— Tanh 



^ ■ ^2Tan-if-)-Tan-i('^-tI-l^V-sgn(2^)K2- (B3) 



For the 2x2 Hamiltonian axx is given by N times 
of Eq. dHIl with Eq. (|B2|. This is consistent with the 
result of Ref. [13! where TV = 2 and E = /i = case IS 
discussed. 

We omit analytical result of the Hall conductivity a^y 
for the bilayer system, because of its lengthincss. Eor 
axy of the 2x2 Hamiltonian ([3]) , the following two terms 
should be added in tr (• • • ) of the polarization function 
(|15p . since 7^ has momentum dependence, 



GlxG+-fyG+idk^-fy) + QlxQ+{dk^ly)Gly 



(B4) 



In this case, finite cutoff is needed to make the Uxy finite. 
Moreover, the Drude-Zener-like behavior for the high- 
energy region cannot be described. As a result, it turns 
out that Hamiltonian ^ is not appropriate to discuss 
the Hall conductivity in a weak-magnetic field. 

The magnetic susceptibility (|18p for bilayer system 
with 4x4 Hamiltonian is obtained as 



127r2c2 



dxf{x)B2{x) 



(B5) 



and 



B2{x) =Im 



{x + iE)2 - i2 



t{x + lE) 



log 



t + x + \V 



-t + x + \T 



(B6) 



In the limit E — > 0, Eq. (|B5p becomes the result obtained 
by Safrau)^ 



12ttcH 



fi-t) - fit) + W 1' dxl^""^ 



(B7) 



0(t-|M|)(l + 31og|^|) (T = 0), (B8) 



127rc2t 



where V means the Cauchy's principal value. On the 
other hand, in the limit t — > 0, we obtain a value two 
times of that of the monolayer system (fTQ]). 

The derivation of the formula (|18p for the 2x2 Hamil- 
tonian ([3]) can be done as in the same way of Ref. In 
this calculation, we should extract terms of the thermo- 
dynamic potential which is proportional to squire of the 
magnetic field in the following representation: 



^2 ^ 2g+g- 



-A+A^ - q\A^_ - q^_A\ 



(B9) 



The formula for N = 2 was also derived in Ref. |30| using 
1x1 ly which has more complicated expression. These 
two expressions can be shown to be equivalent using some 
identities. To calculate susceptibility for Eq. (I3|), a finite 
cutoff is needed to make the susceptibility finite. This 
result shows diamagnetism near the zero energy point, 
but it is not appropriate for quantitative analysis. 
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